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In  this  paper  we  present  a  conceptual  design  and  parameters  for  a  compact, 
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wavelength  to  be  2.141  pm  because  of  propagation  and  generation  advantages  as  well  as 
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MW  and  the  input  is  provided  by  a  low  average  power  (<  50W)  FEL  oscillator. 
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grazing  angle  relay  mirror  to  be  close  to  the  exit  of  the  wiggler  without  exceeding  the 
mirror  damage  intensity  threshold  level,  which  is  taken  to  be  50  kW/cm^.  In  this  design 
the  wiggler  length  is  <  2  m  and,  depending  on  the  mirror  grazing  angle  and  the  optical 
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L  Introduction 

The  free  eleetron  laser  (PEL)  is  capable  of  producing  high  average  power  at  high 
efficiency  without  the  conventional  thermal  management  and  waste  issues  associated 
with  other  laser  systems  [1].  In  addition,  the  operating  wavelength  can  be  chosen  for 
optimized  propagation  in  a  maritime  environment  [2].  These  unique  features  make  the 
PEL  a  leading  candidate  for  naval  directed  energy  applications  [3-6], 

Propagation  of  high  power  laser  beams  in  the  atmosphere  involves  a  complex, 
linear  and  nonlinear  interaction  between  various  physical  processes.  The 
theoretical/numerical  model  used  in  this  study  includes  the  effects  of  aerosol  and 
molecular  scattering,  aerosol  heating  and  vaporization,  thermal  blooming  due  to  aerosol 
and  molecular  absorption,  atmospheric  turbulence,  and  beam  quality.  These  processes  are 
modeled  in  a  fully  three-dimensional  and  time-dependent  manner.  It  is  found  that 
aerosols,  which  consist  of  water,  sea  salt,  organic  matter,  dust,  soot,  biomass  smoke, 
urban  pollutants,  etc.,  are  particularly  important  because  they  result  in  laser  scattering, 
absorption  and  enhanced  thermal  blooming. 

In  this  paper  we  present  parameters  for  a  conceptual  design  of  a  compact, 
megawatt-class,  optically  guided,  free-electron  laser  amplifier.  There  are  a  number 
important  characteristics  and  advantages  of  the  proposed  configuration.  In  an  optically 
guided  EEL  amplifier  the  filling  factor,  i.e.,  the  ratio  of  the  electron  beam  to  the  optical 
beam  cross-sectional  area,  remains  relatively  constant  throughout  the  wiggler  region,  and 
the  radiation  amplitude  increases  exponentially  until  saturation.  The  optical  beam  upon 
exiting  the  wiggler  can  be  pinched  by  focusing  the  electron  beam  with  external  focusing 
fields.  The  Rayleigh  range  of  the  pinched  optical  beam  can  be  short.  As  a  result  the  first 
relay  mirror  can  be  placed  close  to  the  wiggler  without  exceeding  the  intensity  damage 
threshold.  The  first  relay  mirror  is  a  grazing  incidence  reflector.  The  EEL  efficiency  can 
be  increased  by  tapering  the  wiggler.  Alternatively,  in  a  uniform  wiggler,  enhanced 
efficiency  can  be  obtained  by  frequency  detuning  the  input  signal.  Eor  the  power  levels 
of  interest  we  choose  the  EEL  wavelength  to  be  2.141  pm  because  of  propagation  and 
generation  advantages  as  well  as  for  eye  safety  reasons. 

Recently  a  high-gain  EEL  amplifier  operating  in  the  optically  guided  regime  has 
been  experimentally  demonstrated  at  the  Brookhaven  National  Laboratory  [7].  Our 


2 


analysis  of  this  experiment  is  in  close  agreement  with  the  experimental  results  and  will  be 
discussed  in  more  detail  in  Sec.  IV. 

II.  Laser  Propagation  in  a  Maritime  Environment 

The  optimum  wavelength  for  efficient  HEL  propagation  depends  on  the 
atmospheric  conditions  and  a  number  of  inter-related  physical  processes  which  include: 
thermal  blooming  due  to  aerosol  and  molecular  absorption  [8],  turbulence  [9],  aerosol 
and  molecular  scattering  [10],  thermal  scattering  due  to  heated  aerosols,  and  aerosol 
heating  and  vaporization  [11-14].  The  relative  importance  of  these  processes  depends  on 
the  parameters  of  the  atmospheric  environment  which  can  vary  significantly  depending 
on  location  and  time. 

Atmospheric  environments  contain  various  types  and  concentrations  of  aerosols. 
Aerosols  can  absorb  laser  energy  and  conductively  heat  the  surrounding  air,  resulting  in 
an  increase  in  thermal  blooming  of  the  HEL  beam  [2,15].  In  general,  aerosols  consist  of 
hygroscopic  and  non-hygroscopic  particles  of  various  sizes  and  chemical  compositions. 
Oceanic  aerosols  consist  of  sea  salt,  water,  and  organic  material.  Non-hygroscopic 
aerosols  are  composed  of  dust,  soot,  biomass  smoke,  and  other  carbon-based  compounds. 
These  aerosols  typically  have  much  larger  absorption  coefficients  than  water-based 
aerosols.  While  they  are  normally  present  in  continental,  rural  and  urban  environments, 
dust  aerosols  can  also  be  present  in  maritime  environments  hundreds  of  miles  from  shore 
[16]. 

i)  Maritime  Atmosphere  Model 

The  Advanced  Navy  Aerosol  Model  (ANAM)  is  used  to  model  the  near-surface 
maritime  environment  [17].  The  ANAM  aerosol  distribution  is  comprised  of  various 
modes  which  represent  aerosols  of  different  compositions  and  sizes.  These  aerosol 
modes  will  absorb  laser  energy  and  vaporize  at  different  rates.  The  aerosol  size 

4 

distribution  function,  F(R,t)  =  '^Fj(R,t) ,  where  R  is  the  aerosol  radius,  is  represented 

;  =0 

as  a  superposition  of  five  “modes”  with  each  mode  representing  aerosols  with  a  particular 
physical  composition  and  origin.  The  aerosol  distribution  can  evolve  with  time,  due  to 
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vaporization  for  example,  the  aerosol  absorption  and  scattering  coefficients  are  also  time- 
dependent.  Mode  0  represents  dust  particles  of  continental  origin,  mode  1  represents 
water-soluble  aerosols,  and  modes  2-4  represent  marine  aerosols  (sea  salt  and  water)  that 
result  from  different  processes. 

We  use  MODTRAN  to  calculate  the  molecular  absorption  and  scattering 
coefficients  [18].  Figure  1  plots  the  gross  scattering  and  absorption  coefficients  due  to 
aerosol  and  molecular  contributions  in  the  water  vapor  transmission  windows, 

Jl  =  1.045, 1.625  and  2.141)j.m .  The  gross  scattering  and  absorption  coefficients  that  we 
obtain  are  comparable  with  in  situ  measurements  [19,20]. 


ii)  HEL  Propagation  Simulations 

To  simulate  HEL  propagation  we  use  the  High  Energy  Laser  Code  for 
Atmospheric  Propagation,  HELCAP  [2,21],  developed  at  the  Naval  Research  Laboratory. 
HELCAP  includes  the  effects  of  i)  aerosol  and  molecular  scattering,  ii)  aerosol  heating 
and  vaporization,  iii)  thermal  blooming  due  to  both  aerosol  and  molecular  absorption,  iv) 
atmospheric  turbulence,  and  v)  laser  beam  quality,  i.e.,  all  of  the  effects  discussed  in  the 
preceeding  section.  It  is  the  first  HEL  propagation  model  which  integrates  all  these 
physical  processes  in  a  fully  three-dimensional,  time-dependent  manner.  In  modeling  the 
aerosol  effects,  we  account  for  the  aerosol  distribution  and  the  various  aerosol  modes 
(water-based,  dust,  soot,  etc.). 

In  HELCAP  the  laser  electric  field  is  represented  as 
Yj  =  A{x,y,z,t)e.x^[  i{coz ! c  -  6>/)]  /2  +  c.c.,  where  (y  =  2;rc//l  is  the  laser 

frequency,  is  a  unit  polarization  vector  in  the  x  direction,  A(x,y,z,t)  is  the  complex 

laser  amplitude  and  the  laser  intensity  is  I  =  c  A  A*  l?>n  .  HELCAP  solves  a  nonlinear 
Schrodinger-like  equation  which  has  the  form 


dA 
d  z 


ic 

2  CO 


VIA 


i  —  {Sn^  +  Sn^g ) 
c 


+  /?) 


A, 


where  a  =  \s  the  total  absorption  coefficient,  P  =  is  the  total 

scattering  coefficient,  and  Sn.,.  and  dn-jg  denote  the  refractive  index  variation  due  to 
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atmospheric  turbulence  and  thermal  blooming  respectively.  )  is  the  molecular 

(aerosol)  absorption  coefficient.  )  is  the  molecular  (aerosol)  scattering 

coefficient.  The  quantities  5n-^g  ,  a  ,  and  p  are  space  and  time-dependent  and 
determined  self-consistently  in  the  presence  of  aerosol  heating  and  vaporization. 

The  transmitted  power  at  the  source  is  denoted  by  P  and  the  laser  is  focused  onto 
a  remote  target  at  a  range  of  5  km  (see  Fig.  2).  The  target  is  taken  to  be  circular  with  an 
area  of  1 00  cm'.  The  propagation  direction  is  along  the  z-axis  and  a  uniform  transverse 
wind,  with  velocity  =  5  m/sec  is  directed  along  the  y-axis.  Atmospheric  turbulence  is 

modeled  by  a  Kolmogorov  spectrum  with  structure  constant  Cl  =10''^  m  '^^ .  The 

pointing  jitter  associated  with  the  laser  beam  is  taken  to  have  an  angular  spread  of  2  prad 
and  a  white  noise  temporal  spectrum. 

Figure  3  plots  the  average  power  on  target  versus  the  transmitted  power,  P ,  for 
the  three  wavelengths  of  interest.  Our  results  show  that  for  a  maritime  environment,  the 
optimum  wavelength  depends  on  the  transmitted  power.  For  P  <  1.5  MW,  propagation 

is  mostly  affected  by  aerosol  scattering  and  the  average  power  on  target  increases  with 
P  .  In  this  regime,  the  1.625  pm,  and  2.141  pm  wavelengths  provide  slightly  greater 
power  on  target  than  1 .045  p  m.  This  is  due  to  the  lower  aerosol  scattering  coefficient 
associated  with  the  longer  wavelengths.  For  P  <  1 .5  MW,  the  propagation  efficiency  is 
roughly  50%  for  the  three  wavelengths  considered.  For  example,  P  =  1  MW  results  in 

(/Jar^e.)  ~  0-55  MW  for  A  =  1.625pm  and  T  =  2.141pm.  However,  for  P  >  1.5  MW, 

thermal  blooming  becomes  important.  In  this  high  power  regime  the  optimum 
wavelength  is  1 .045  pm  due  to  the  lower  molecular  absorption  coefficient  in  that  water 

vapor  window.  For  P  =  3  MW,  ^FJarsei)  ~  ^  f  045  pm  while  ^P,arget)  ~  0.8  MW 

for  1.625  pm  and  2.141  pm. 

For  the  power  levels  of  interest,  i.e.,  P  <  1 .5  MW  average  power,  the  FEL 
wavelength  is  chosen  to  be  2.141pm .  For  maritime  propagation  there  is  no  distinct 
advantage  in  operating  the  FEL  at  either  1.045  pm  or  1 .625  pm .  In  the  following 
Section,  we  present  design  parameters  for  a  high-gain  FEL  amplifier  operating  at 
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2.141)xm  capable  of  achieving  ~  1.5  MW  of  average  output  power.  We  note  that  for  the 
intensity  levels  achievable  in  present  HEL  systems,  it  is  the  average  power,  and  not  the 
peak  power,  that  determines  the  propagation  and  lethality  characteristics  of  the  HEL 
beam. 


III.  Optically  Guided  FEL  Amplifier  Operating  at  2.141  pm 

Using  the  Source  Dependent  Expansion  (SDE)  formulation  [22],  expressions  for 
the  optical  growth  rate,  wavenumber  shift,  spot  size  and  optical  curvature  of  the  optical 
beam  in  an  FEL  amplifier  are  obtained  in  the  Appendix.  These  results  are  used  to  design 
a  compact  MW  FEL  in  the  optically  guided,  high  gain  regime,  i.e.,  having  a  matched 
optical  spot  size.  In  addition,  focusing  the  electron  beam  at  the  wiggler  exit  is  shown  to 
result  in  a  pinched  optical  beam  with  a  shortened  Rayleigh  range.  A  schematic  diagram 
of  the  FEL  amplifier  configuration  is  shown  in  Fig.  4.  Finally,  it  is  shown  that  by 
frequency  detuning  the  FEL  efficiency  can  be  more  than  doubled.  The  following 
discussion  is  based  on  the  design  parameters  listed  in  Tables  I  -  II  and  for  the  case  of  no 
frequency  detuning. 

i)  Operating  Wavelength 

In  the  present  design  the  resonant  operating  wavelength  of  a  linearly  polarized 
optical  beam  is  A  =  (1  +  Ai W2)/2y"  =  2.141pm,  where  y  =  T'^[MeV]/0.51 1  +  1  is 

the  total  relativistic  mass  factor  (see  Tables  I  and  II). 


ii)  Gain  Length 

The  power  gain  length  in  the  optically  guided  regime,  in  the  absence  of  frequency 
detuning  (Aru  =  0) ,  is 


'  (2  +  3/) 

y 

1/2 

K 

2(1  +  2/)" 

V 

z 

V 

^  n  ^ 

=  14cm 


(1) 


where  /  =  /  Rj  <1  is  the  filling  factor,  v  =  /4[kA]/17  is  Budker’s  parameter,  4  is 

the  electron  beam  current  in  kA,  K  =  q  Bfy  /iTrmc^  =1.9  is  the  wiggler  strength 
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parameter,  =  5kG  is  the  wiggler  field  and  =  4  cm  is  the  wiggler  period.  The 
transverse  profiles  of  the  optical  and  electron  beam  are  Gaussian  with  radii  Rf  and  Rj, , 
respectively  [see  Eqs.  (A2)  and  (A4)]. 


iii)  FEL  Efficiency 


In  the  optically  guided,  high  gain  regime  the  efficiency  of  conversion  from 
electron  beam  power  to  optical  power  is  given  by  [3,23] 


ti  = 


i.M  _  /‘'+2  +  3/)”^ 

V 

1/2 

2^  (1  +  2/)' 

K 


+  K"  n 


=  0.8%. 


(2) 


Here,  AA:  is  the  wavenumber  shift,  defined  in  Eq.  (A3).  Free-electron  laser  efficiency 
can  be  increased  by  tapering  the  wiggler  [4,  23,  24]  or  by  frequency  detuning  [3,23]. 


iv)  Optical  Guiding 

The  FEL  optical  beam  can  undergo  refractive  guiding  within  the  wiggler.  An 
envelope  equation  for  the  optical  beam  spot  size  can  be  obtained  by  using  the  results  in 
the  Appendix,  and  the  condition  for  a  matched  optical  beam  can  be  obtained.  For  optical 
guiding  the  power  gain  length  is  nearly  equal  to  the  free  space  Rayleigh  range,  that  is 

,  (2  +  3/)"' 


where  =  n  R^  !  X  =  22  cm  is  the  Rayleigh  range.  The  optical  guiding  condition  in  Eq. 
(3)  can  be  rewritten  in  terms  of  the  beam  current, 

.-2  /^(l  +  2/)^ 


I,[kA\  =  l.lxlO' 


(2  +  3/)' 


'  (l  +  K-/2y 

^  y 

=  IkA. 


(4) 


In  addition,  the  centroid  of  the  FEL  optical  beam  can  be  steered  by  gradually  bending  the 
electron  beam. 
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v)  Saturated  Peak  and  Average  Optical  Power 

The  peak  optical  power  at  saturation  is  /*  =  77  =1  GW  ,  where  =  81  MV 

is  the  electron  beam  voltage.  The  average  power  at  saturation  is 


(P)=r]DIJ,  =  1.5MW, 


(5) 


where  D  =  1 .5 x  1 0  ^  is  the  duty  factor  and  r]  =  \.1%  ( hco! o  =  -1.5%). 

vi)  Electron  Pulse  Slippage  (Lethargy) 

In  the  high-gain  regime  electron  slippage  (lethargy)  is  shown  to  be  significantly 
less  than  in  the  low  gain  regime  and  is  not  a  concern  for  this  design.  The  slippage  length, 
i.e.,  the  separation  between  the  optical  and  electron  pulse,  is  5  s  (v^  -  v.)  / c ,  where 

v^  =  daldk  is  the  optical  pulse  group  velocity,  v,  =  c{\  -  Mlyl)  is  the  axial  electron 

beam  velocity,  and  y,  =  y  I -Jl  -1-  {K^  /  2)  .  Making  use  of  the  high  gain  FEL  dispersion 

relation  [3,23],  the  group  velocity  is  v^  s  c  (l  -  1  / 3  y?  -{X! n  R,  y  /  2),  and  the  slippage 


length  is 


where  I is  the  number  of  wiggler  periods.  The  well  known  slippage  length 

in  the  low  gain  regime  is  X  N ^  .  Slippage  in  the  high  gain  regime  is  reduced  by  more 
than  a  factor  of  three. 

vii)  Pinched  Optical  Beam 

The  intensity  on  the  first  relay  mirror  must  be  below  a  damage  threshold  level. 

To  further  reduce  the  intensity  on  the  mirror  the  optical  beam  can  be  pinched  at  the 
wiggler  exit  in  order  to  shorten  the  Rayleigh  range  [3,4].  By  employing  external  focusing 
fields  near  the  wiggler  exit  the  electron  beam  and  in  turn  the  optical  beam  can  be 
pinched.  Making  use  of  grazing  incidence  geometry  the  intensity  on  the  mirror  can  be 
further  reduced. 
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The  spot  size  of  the  optical  beam  (z)  [obtained  by  solving  Eqs.  (A5),  (A6)  and 
(A8)]  and  the  radius  of  the  electron  beam  7?^  (z)  are  plotted  as  functions  of  the 
propagation  distance  z  in  Fig.  5.  The  entrance  to  the  wiggler  is  atz  =  0  and  extends 
to  z  =  =  1 .5  m  .  In  Fig.  5(a)  the  radius  of  the  electron  beam  is  constant  and  the  optical 

beam  is  guided  through  the  wiggler  with  a  constant  spot  size.  Upon  exiting  the  wiggler, 
the  optical  beam  diffracts.  In  the  example  of  Fig.  5(b)  the  electron  beam  is  focused  by 
external  fields  applied  near  the  wiggler  exit.  This  pinches  the  optical  beam  and  as  a 
result,  the  optical  beam  diffracts  more  rapidly  than  in  the  unpinched  example. 


viii)  Optical  Beam  Quality  and  Beam  Spreading 

It  is  common  practice  to  characterize  the  higher  order  modal  content  of  a  laser 
beam  by  a  beam  quality  parameter  denoted  by  MT  The  quantity  >  1  is  a  “times 
diffraction-limited”  parameter  which,  for  a  fundamental  Gaussian  beam,  is  unity.  Upon 
exiting  the  FEL  interaction  region,  i.e.,  wiggler,  the  spatial  spreading  of  the  optical  beam 
is  given  by 

W\z)  =  l{r,z)2Krdr  =  If/  +  el{z  -  (7) 

^  0 


where  W  is  the  laser  beam  spot  size,  P  is  the  power  and  1  is  the  intensity.  In  Eq.(  7),  fr), 
is  the  minimum  spot  size  at  z  =  =  X/{n:W^)  is  the  diffraction  (spreading)  angle 

of  an  ideal  Gaussian  beam  and  MMs  a  measure  of  the  beam  quality.  For  a  non-ideal 

optical  beam  the  diffraction  angle  is  6  at  the  exit  of  the  wiggler.  Using  the 
paraxial  wave  equation  in  Sec.  II,  it  can  be  shown  that  the  general  expression  for  the 
evolution  of  is 


M^(z)  = 


An" 


UU  \ 

[ir^dr  I 


0 


5/ 


1/2 


dr 


dr 


r  dr 


>  dr 


1/2 


(8) 
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where  the  complex  laser  amplitude  has  been  expressed  in  terms  of  the  intensity,  /,  and 
phase,  ^,as  A  ^  (r,z)exp(i(/){r,z)) .  When  the  laser  beam  propagates  in  a  region 

of  uniform  refractive  index,  i.e.,  outside  of  the  wiggler,  M^is  constant.  Simulations 
using  the  PEL  code  MEDUSA  indicate  that  the  beam  quality  of  an  optically  guided 
PEL  beam  is  fairly  good.  Por  parameters  similar  to  those  in  Table  I  and  II ,  simulation 
results  indicate  the  «  1 .5  at  the  exit  of  the  wiggler  [25] 


ix)  Distance  Between  First  Relay  Mirror  and  Wiggler 

To  avoid  damage,  the  distance  between  the  wiggler  and  relay  mirror  must  exceed 


L 


Z. 


relay 


,1/2 


sin''"  0 


(9) 


^  \  damage  j  j 
2 


where  0  is  the  tilt  angle  (Pig.  4),  -  n  I Z  \s.  the  Rayleigh  range  at  the  exit  of  the 

wiggler,  (/)  =  2(P)/  t:  Rl  is  the  average  intensity  at  exit  of  the  wiggler,  and 


(-^damage)  ~  kW/cm ^  is  the  average  intensity  damage  threshold  of  the  mliTor.  Using 
the  parameters  in  Tables  I,  II,  and  III,  we  obtain  for  an  unpinched  optical  beam. 


L 


relay 


and  for  a  pinched  beam 


L 


relay 


15m,  0  =  45° 

8m,  0  =  5° 

5  m,  0  =  45° 

3m,  0  =  5° 


(10) 


(11) 


x)  Frequency  Detuning  and  Efficiency  Enhancement 

In  an  PEL  the  efficiency  can  be  increased  by  operating  at  a  frequency  that  is 
detuned  from  the  resonant  value  cOp_~2ncl  Z .  In  general,  however,  operating  off- 
resonance  leads  to  a  reduction  in  the  growth  rate  [1].  If  the  operating  frequency  is 
denoted  by  co ,  the  modification  to  the  present  analysis  is  the  replacement  of  the 
wavenumber  shift  AA:  with  AA  -  (fi;-(5),^)/(2cy_?)  in  the  expression  for  the  efficiency  in 
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Eq.  (2)  and  in  Eq.  (A7)  which  leads  to  the  dispersion  relation.  Note  that  the  left  hand  side 
of  Eq.  (A8)  is  not  changed. 

IV.  Discussion  and  Conclusions 

Conceptual  design  parameters  for  a  compact  MW-class  FEE  amplifier  operating 
in  the  optically  guided,  high  gain  regime  are  presented.  The  design  presented  here 
assumes  that  the  FEE  amplifier  is  driven  by  a  high  quality  RE  linac  generated  electron 
beam.  Optical  guiding  allows  the  interaction  length  to  be  sufficiently  long  to  achieve 
substantial  growth  of  the  input  signal.  Electron  pulse  slippage,  which  can  limit  the  FEE 
interaction  length,  is  shown  to  be  significantly  reduced  for  a  high  gain  amplifier.  By 
pinching  the  optical  beam  at  the  wiggler  exit  the  Rayleigh  range  is  shortened,  thus 
reducing  the  required  distance  to  the  first  relay  mirror.  The  optical  beam  can  be  pinched 
by  focusing  the  electron  beam  with  external  magnetic  fields  or  by  utilizing  the  betatron 
oscillation  of  the  electron  beam  envelope.  Grazing  incidence  geometry  permits  a  further 
reduction  in  intensity  to  avoid  damage  to  the  mirror. 

In  this  design  the  operating  wavelength  is  2.141  pm  and  is  chosen  because  of 
propagation  and  generation  advantages  as  well  as  for  eye  safety  reasons.  For  a  1  kA,  81 
MV  electron  beam,  a  wiggler  field  of  5  kG,  and  a  frequency  detuning  of  - 1 .5%  ,  the 
power  gain  length  is  14  cm,  the  intrinsic  efficiency  is  ~  1.2%,  the  optical  output  power  is 
1  GW,  and  the  average  output  power  is  1.5  MW,  assuming  a  duty  factor  of  1.5x  10“^  for 
the  RF  linac.  Focusing  the  electron  beam  radius  down  by  a  factor  ~  3  at  the  exit  of  the  ~ 

1 .5  m  wiggler  reduces  the  optical  beam  Rayleigh  range  by  a  factor  of  ~  4.  Depending  on 
the  mirror  grazing  angle  and  the  optical  beam  quality,  the  distance  from  the  wiggler  exit 
to  the  first  relay  mirror  can  be  <  3  m.  The  length  of  the  optical  system,  which  includes 
the  amplifier  section  and  the  first  relay  mirror,  can  be  <  5  m. 

A  recent  experiment  at  the  Brookhaven  National  Eaboratory  (BNE)  has 
demonstrated  the  optically  guided,  high-gain  FEE  amplifier  concept.  This  experiment 
utilized  a  10  m  long  wiggler  and  an  S-band  RF  linac  that  produced  a  100  MeV,  3 00 A 
electron  pulse  of  radius  0.2  mm.  The  wiggler  period  and  strength  parameter  were  3.9  cm 
and  1.1,  respectively.  Experiments  showed  optical  guiding  in  a  high  gain  FEE  amplifier 
operating  at  T  =  0.8  //  w  using  a  Ti-Sapphire  laser  as  the  seed.  Dipole  steering  magnets 
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located  along  the  wiggler  were  used  to  deflect  the  electron  beam  by  ~  Imrad  and  the  light 
beam  was  shown  to  be  steered  and  guided  by  the  electron  beam.  The  transverse  profile  of 
the  optically  guided  light  beam  was  shown  to  be  less  than  half  the  size  of  an  unguided 
diffracting  light  beam.  The  power  peak  output  power  was  ~150  MW  [7].  Our  simulations 
of  the  BNL  experiment  predict  a  power  e-folding  length  of  42  cm,  a  guided  spot  size  of 
0.24  mm  (1/e  of  intensity),  and  a  peak  power  output  of  ~180  MW,  in  close  agreement 
with  the  experimental  results. 

Simulations  of  the  PEL  amplifier,  taking  into  account  emittance,  energy  spread, 
and  non-uniformity  of  the  optical  guiding  due  to  variation  of  the  electron  current  within 
the  pulse,  are  required  to  fully  evaluate  the  present  design.  Also,  while  estimates  indicate 
that  a  modest  external  magnetic  field  is  sufficient  to  focus  the  beam,  more  refined 
analysis  is  required  to  design  the  pinch  region  of  the  wiggler.  In  the  absence  of  optical 
pinching,  the  distance  to  the  first  relay  mirror  can  be  <  10  m. 
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Appendix 


The  source  dependent  expansion  (SDE)  formulation  [22]  can  be  used  to  analyze 
the  dynamics  of  the  optical  beam  in  FELs.  A  summary  of  the  principal  SDE  equations  is 
given  in  this  appendix.  For  a  planar  wiggler  the  electric  field  of  the  radiation  can  be 
taken  to  be  linearly  polarized,  and  given  by 

E{r,z)  =j  A{r,z)  exp{ico{z I c  -  +  c.c.  (Al) 

For  the  fundamental  Gaussian  mode  the  complex  envelope  E„(r,z)  is  written  as 

A(r,z)  =  A„exp{iA^(z))exp{-(\  -  ia{z))r^  /  RI(z)),  (A2) 

where  A^^  is  the  amplitude  on  axis, 


T^(z)  =  [  (M(z')-  /r(z'))^',  (A3) 

Ak(z)  is  the  wavenumber  shift,  r(z)  is  the  spatial  growth  rate,  a(z)  is  related  to  the 
curvature  of  the  optical  wavefronts  and  (z)  is  the  spot  size  of  the  beam.  To  allow  for 
an  electron  beam  with  a  radius  R/,  (z)  that  varies  along  the  propagation  axis,  the  electron 
beam  density  is  taken  to  be 

n,iz,r)  =  n,(0)(Rli0)/ RUz)y^p(-r^^  / RU^))-  (A4) 

The  SDE  equations  for  the  spot  size  and  a(z)  are  given  by  [3] 
dR,(z)  2ca(z)  1 


and 


d  z 


da{z) 


CO 


=  -H,  {z)R,Xz), 


d  z 

respectively.  Here, 


2c)l  +  «  (z)  ^ 

- ,  =  2(//^(z)  -  a{z)H,{z)), 

^co)  Rliz) 


H{z)  =  4 


V  c 


y  CO 


2  /,2  ^ 


K^k 


\  +  12 


f(^) 


Rliz)  (l  +  2f{z)Y  {Ak(z)  -  /r(z))- 


(A5) 


(A6) 


(A7) 
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where  u  =  4[lcA]/17  is Budker’s parameter, =  IkIX^,  f{z)  -  RI{z)I R\{z)  is 


the  filling  factor  and  the  subscripts  R  and  I  denote  real  and  imaginary  parts.  The  final 
SDE  equation  for  the  phase  shift  and  growth  rate  is  given  by 


A^(z)  -  /r(z)  + 


c  1  (1  +  a-{z))  .{\-ia{z))  dRi(z)  ^  I  dajz)  ^  ^  f(z))H(-) 


yco)  Rliz) 


R,  (z)  d  z 


(A8) 

Equations  (A5),  (A6)  and  (A8)  describe  the  optical  beam  in  the  high  gain  EEL  amplifier. 
The  electron  beam  can  be  focused  to  produce  a  pinched  optical  beam. 
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Figure  1:  Aerosol  scattering,  aerosol  absorption,  and  water  vapor  (molecular)  absorption 
coefficients  used  for  the  model  maritime  environment  of  the  propagation  simulation. 
Aerosol  scattering  and  absorption  are  calculated  using  ANAM,  water  vapor  absorption  is 
calculated  using  MODTRAN4. 
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Figure  2;  Schematic  of  laser  and  target  configuration  used  in  the  full 
scale  propagation  simulations. 
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Figure  3:  Average  power  on  target,  (/’orget )  ?  versus  transmitted  power,  P  ,  in  a 
maritime  environment  for  the  wavelengths  X  =  1.045, 1.625,  2.141  pm  .  Initial  beam 
profile  has  -  50  cm,  £)q  =  80  cm .  Simulation  geometry  is  shown  in  Fig.  2.  Target 
range,  L  =  5  km;  beam  focus  =  5  km;  target  area  =100  cm  ;  wind  speed,  =  5 
m/sec;  turbulence  strength,  =  10“'^  m'"^ ;  pointing  jitter  angular  amplitude  =  2 
prad  (white  noise). 
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Figure  4.  Schematic  of  high  gain  FEL  amplifier  with  a  grazing  relay  mirror.  The  input 
signal  can  be  obtained  a  low  average  power  FEL  oscillator.  The  radiation  beam  is 
optically  guided  in  the  wiggler  and  optically  pinched  at  the  exit.  The  pinched  optical 
beam  has  a  shortened  Rayleigh  range  and  undergoes  rapid  diffraction  upon  exiting  the 
wiggler.  Employing  a  grazing  incidence  configuration  the  resultant  footprint  on  the  relay 
mirror  can  be  made  sufficiently  large  to  avoid  damage. 
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Figure  5.  Plots  of  optical  beam  radius  (Ri)  and  electron  beam  radius  (  i?,, )  versus 
propagation  distance  z  .  Wiggler  entrance  is  at  z  =  0  and  the  exit  is  at  z  =  ~  1 .5  m  . 

The  FEL  operates  in  the  high-gain  regime,  where  the  optical  beam  maintains  a  constant 
spot  size  by  optical  guiding.  Outside  the  wiggler  the  optical  beam  diffracts  and  spreads 
transversely,  as  indicated  in  (a).  By  focusing  the  electron  beam  at  the  wiggler  exit  the 
optical  beam  can  be  pinched  down  to  pinch  =  ^.2 1  mm  as  indicated  in  (b).  The 
pinched  optical  beam  diffracts  more  rapidly  compared  to  the  unpinched  example,  i.e.,  has 
a  shorter  Rayleigh  range,  ^  R^pinch  =  6.7  cm .  In  generating  this  figure  a 

frequency  detuning  of  Aco/  co  -  -1.5%  was  used. 
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Table  I.  Wiggler  Parameters 


Wiggler  Field, 

5kG 

Wiggler  Parameter,  K 

1.86 

Period, 

4  cm 

Length, 

1.5  m 

Table  II.  Electron  Beam  Parameters 


Energy, 

81  MeV 

Current, 

1  kA 

Radius,  7?^ 

0.02  cm 

Normalized  Emittance, 

15  mm-mrad 

Duty  Factor,  D 

1.5  xlO”' 

RF  Linac  Frequency 

750  MHz 

Micro  Pulse  Duration 

2  psec 

Betatron  Wavelength, 

X,=2yXJK 

7m 
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Table  III.  FEL  Optical  Beam  Parameters 


Operating  Wavelength,  X 

2.141  pm 

2.141  pm 

Frequency  Detuning,  ^.colco 

0 

-1.5% 

Spot  Size  in  Wiggler,  R, 

0.33  mm 

0.4  mm 

Filling  Factor,  f 

0.36 

0.25 

Rayleigh  Range  (unpinched,  at  wiggler  exit ), 

16  cm 

23  cm 

Gain  Length  (power), 

14  cm 

14  cm 

Efficiency  (uniform  wiggler),  77 

O.ITA 

1.2% 

Average  Intensity  (unpinched,  at  wiggler  exit ),  (/) 

0.53  GWW 

0.58  GWW 

Pinched  Laser  Spot  Size, 

0.17  mm 

0.21  mm 

Rayleigh  Range  (pinched,  at  wiggler  exit), 

4  cm 

6.7  cm 

Average  Pinched  Intensity,  ^ 

2.2  GWW 

2.1  GW/cm^ 

Peak  Output  Power,  P 

0.63  GW 

1  GW 

Average  Output  Power,  [p'^ 

1  MW 

1.5  MW 

Number  of  Power  e  -  Folds, 

10.6 

10.5 

Average  Input  Power,  ) 

43  W 

40  W 

Optical  Beam  Quality, 

1.5 

1.5 
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Table  IV.  First  Relay  Mirror  Parameters  (for  Aco!  co  =  -1.5%) 


Tilt  Angle,  0  [degrees] 

45° 

5° 

Distance  from  Wiggler  (unpinched),  [m] 

14 

5 

Distance  from  Wiggler  (pinched),  [m] 

8 

3 

The  average  intensity  damage  threshold  was  taken  to  be 
{l^.)  =  50kW/cm= 
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